Synchronization of Coupled Anizochronous Auto-Oscillating Systems 
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The particular properties of synchronization are discussed for coupled auto-oscillating systems, 
which are characterized by non-quadratic law of potential dependence on the coordinate. In par- 
ticular, structure of the parameter plane (frequency mismatch - coupling value) is considered for 
coupled van der Pol - Duffing oscillators. The arrangement of synchronization tongues and the 
particular properties of their internal structure in the parameter space are revealed. The features 
of attractors in the phase space are discussed. 
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I. INTRODUCTION 

The system of two coupled van der Pol oscillators is 
one of the basic models of nonlinear dynamics demonstra- 
ting the phenomenon of mutual synchronization. There 
are many papers on this theme [see 1-15] because 
this system demonstrates a lot of interesting oscilla- 
tion regimes and types of behavior, such as synchronous 
and quasiperiodic regimes, synchronization with different 
phase conditions [1, 15], the "oscillation death" effect [1, 
2], global bifurcations [7], etc. Moreover, the results of 
this system investigation may be used in the analysis of 
different electronic, biological and chemical systems [3 - 
5, 11 - 13]. Interest in this problem does not decrease 
because of its significance and dynamics variety [see 14]. 

For the van der Pol system is only single non-linearity 
characteristic. It is non-linear dissipation, which limits 
amplitude of auto-oscillation. Potential function is the 
same as in the case of the harmonic oscillator. Deve- 
loping this system, we consider van der Pol - Duffing 
oscillator, for which additional non-linearity of Duffing 
oscillator type is introduced [1, 15]. Potential contains 
in this case a component in the form of the coordinate 
raised to the fourth power. This component takes into 
account an important physical effect, namely an oppor- 
tunity of anisochronous oscillations, i.e. the dependence 
of oscillation frequency on oscillation amplitude [1, 15]. 

One of the most important methods of such systems 
synchronization analysis is the quasiharmonic approxi- 
mation. It consists in the construction of equations for 
the amplitude and the phase, which are slowly changing 
against the background of oscillations with the eigenfre- 
quency of oscillators [1, 15]. But this method is effective 
only when the excess above the threshold of bifurcation 
of auto-oscillation appearance is small. So it designates 
only the structure of the main synchronization tongue. 



At the same time the overall picture of the parameter 
space structure for coupled auto-oscillating systems is 
substantially not yet to be established. (It is in contrast 
to the case of forced synchronization by external har- 
monic signal [1, 16, 17].) It concerns especially the case 
of anisochronous system, when the structure of characte- 
ristic regime areas depends on the additional parameter, 
which is responsible for the form of potential function. 

In the present paper the parameter space structure 
is investigated for coupled auto-oscillating systems with 
several characteristic types of potential function, which 
are leading to the anisochronous oscillations. There are 
such cases considered, when quasiharmonic approxima- 
tion is ineffective or unusable at all. 



II. COUPLED VAN DER POL OSCILLATORS 

The system of differential equations describing the in- 
teraction between auto-oscillating systems is of the form 
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Here A is the parameter characterizing the excess 
above the threshold of Andronov - Hopf bifurcation in 
autonomous oscillators, /i is the coefficient of dissipative 
coupling, Ui(x) and U2{y) are the potential functions of 
each oscillator. Their difference determines non-identity 
of oscillators. 

In simple case of quadratic potential 
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the system (1) is the system of two coupled van der Pol 
oscillators, where parameter 8 is the frequency mismatch 
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between the second and the first oscillators: 

x — (A — x 2 )x + x + fi[i, — y) = 0, 

y - (A - y 2 )y + (1 + % + ^(y - A) = 0. 

To determine the parameter space structure of the 
system (3) we shall use the method of dynamic regime 
chart construction [18]. Within the framework of such a 
method we shall mark the oscillation period of the sys- 
tem of coupled oscillators by means of different colors 
(gray color hues) on the parameter plane (frequency mis- 
match S - coupling value fi) . White color corresponds to 
the chaotic or quasiperiodic motions. Cycle periods were 
calculated by means of the Poincare section method: this 
is the number of points of intersection of the phase tra- 
jectory on the attractor and the surface chosen as the 
Poincare section. Only those crossings were taken into 
account that correspond to the trajectories coming to 
the surface from the one side. The number of points of 
intersection of the phase trajectory and the intersecting 
surface was considered as the oscillation period. 

The system under investigation is characterized by 
four-dimensional phase space (x,x,y,y). Therefore 
three-dimensional hypersurface that is preset by means 
of some additional condition, e.g. zero velocity of the se- 
cond oscillator y = may serve as the Poincare section. 
In that case the number n of points of intersection of tra- 



jectory and the hypersurface was determined. Colors on 
the charts are chosen in accordance with the period n. 

The chart of dynamic regimes obtained in such a way 
for the van der Pol oscillators (3) is given in Fig. Q] (a) 
on the (S, /^)-plane for A = 2.5. The value of the control 
parameter is chosen so large, that traditional quasihar- 
monic approximation [1, 14, 15] is not efficient. 

On the chart in Fig. [T] (a) there are shown: main syn- 
chronization tongue with the period 1; area of quasipe- 
riodic regimes with the embedded system of higher or- 
der synchronization tongues; the area of the "oscillation 
death" effect [1, 2] which corresponds to stabilization of 
the equilibrium state point at the origin due to suffi- 
ciently strong dissipative coupling. 



III. COUPLED VAN DER POL - DUFFING 
OSCILLATORS 

Now we pass on to the van der Pol - Duffing model, 
to which the potentials 

t/iOz) = y+/^p U 2 (y) = (l + S)^+(3^ (4) 
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and dynamic equations in the form of 

x — (A — x 2 )x + x + [3x 3 + fx(x — y) = 0, 

jj-(X- y 2 )y + (1 + S)y + (3y 3 + M (y - x) = 0. (5) 

correspond. 

The parameter j3 defines perturbation quantity of the 
quadratic potential and, correspondingly, is responsible 
for anisochronism of small oscillations [1, 15]. 

On the charts in Fig. Q] (b, c) there are shown the dy- 
namic regime chart for the system (5) and its enlarged 
fragment for j3 = 1 that is the case of sufficiently great 
anisochronism. One can see in Fig. [T] (b) that the pre- 
sence of anisochronism leads to the displacement of syn- 
chronization tongues towards the greater values of fre- 
quency mismatch. Synchronization tongues become so 
wide that one can see the situation of their overlap- 
ping, which is characteristic for standard circle map [18]. 
Internal structure of this tongues changes too: period 
doublings and transition to chaotic dynamics inside them 
are observed. It may be seen clearly in enlarged fragment 
of the dynamic regime chart, shown in Fig. [T] (c). 

With the increase of the parameter of anisochronism (3 
an onset of new " islands" of doubled periods takes place 
inside the synchronization tongues locating at greater 
values of frequency mismatch 8 (see Fig. [2]). Synchroniza- 
tion tongues change their shape with the further increase 
of parameter (3. They fall into two similar in structure 
pieces, inside which structures called "crossroad area" 
exist [18]. These structures are characteristic for systems 
with period doublings. Practically, two systems of syn- 
chronization tongues may be observed in the parameters 
plane: tops of the first one are disposed along the line of 
zero coupling, tops of the second one - along the boun- 
dary of main synchronization area 1, which is the line 
of Neimark - Sacker bifurcation. Chaotic area appears 
between these two systems of synchronization tongues. 

Phase plane portraits of attractors computed at several 
characteristic points within the domains of quasiperiodic 
regimes are given in Fig. [3] on (x, x) and {y, y) planes for 
the first and the second oscillators correspondingly. To 
put it more precisely these portraits are projections of 
the attractor in phase space [x, x, y, y) onto correspon- 
ding planes (x, x) and (y,y). For small values of cou- 
pling parameter it may be seen that portraits of attrac- 
tors look like slightly disturbed limit cycles of individual 
oscillators. At the same time the trajectory never ap- 
proaches the neighborhood of the point of origin. With 
the increase of the coupling parameter the trajectories of 
individual oscillators become more perturbed and there 
exists a kind of a threshold value when the trajectory of 
the first oscillator may achieve the neighborhood of the 
origin. Corresponding phase portrait looks like the area 
completely "filled" with trajectories. Evidently, this si- 
tuation corresponds to the case of ill-determined phase of 
this oscillator. When approaching the oscillation death 
area an essential downsizing of attractors occurs (see 
scales in coordinate axis). The orbit ceases again to 



achieve neighborhood of the origin, but structure of phase 
portraits and shapes of attractors are different from the 
case of small coupling value. 




FIG. 2: Dynamic regime chart for the system (5) for A = 2.5, 

f3 = 3. 



IV. COUPLED OSCILLATORS WITH THE 
POTENTIAL OF THE FOURTH DEGREE 

Now we consider the case when the potential does not 
contain quadratic component and is in the form of 

U 1 (x)=f3 1 ^, U 2 (y)=p 2 ^. (6) 

Corresponding dynamic equations are 

x — (A — x 2 )x + (3\x 3 + fi(x — y) = 0, 
y - (A - y 2 )y + fry 3 + fi(y -x)=Q. 

We shall note that the potential of this type was offered 
by Ueda for common (i.e. not auto-oscillating) oscilla- 
tors. Corresponding model was called Ueda's oscillator 
[18]. 

In this case it is impossible to talk about eigenfrequen- 
cies of oscillators and, therefore, to introduce their mu- 
tual frequency mismatch. But we shall provide for non- 
identity of oscillators by means of different values of the 
parameter [3 1 which is responsible for the shape of the 
potential well (6). If we choose j3\ = 1 for the first oscil- 
lator, we shall vary the parameter (3 2 — (3- 

Structure of the parameter plane {(3, /x) constructed 
for the system (7) is shown in Fig. [3] for A = 1. It is 
remarkable that the system of synchronization tongues 
with different periods exists as before in spite of the fact 
that slightly excited autonomous oscillators have no fixed 
eigenfrequencies of small oscillations at all. Nevertheless, 
different steepness of the potential well of the fourth de- 
gree results in possibility of different rhythms in separate 
oscillators and in synchronization with various periods. 

Absence of the oscillation death area in spite of dissi- 
pative nature of coupling is yet another essentially new 




FIG. 4: Dynamic regime chart for the system (7) for A = 1 and phase plane portraits computed at several chosen points within 
the synchronization tongues with periods 1, 2, 3 and 4. 
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FIG. 5: Dynamic regime chart and phase plane portrait com- 
puted at characteristic point for the system with the potential 
of the fourth degree (10) for A = 2.5. 

moment. This fact may be easily proved by means of 
combined equations (7) . In order to this regime was pos- 
sible, the point of origin x = y — should be attracting. 
System of equations (7) linearized about this point has 
the form of 

x + (n - X)x = fiy, 
y + {H~ X)y = fix. 

Looking for the solution of the system (8) in the form 
of e at , we shall obtain simple characteristic equation 

a 2 (a- X)(a + 2fi- X) = 0. (9) 

One can see that at least one root of this equation 
a = X is always positive when auto-oscillations exist 
and A > 0. So state of equilibrium in the point of ori- 
gin is always unstable. At the same time the condition 
fi = X becomes apparent in the structure of the dynamic 
regime chart. It is the boundary of the main synchroniza- 
tion tongue when j3 3> 1 and oscillators (7) are strongly 
non-identical. One more feature consists in the fact that 
tops of the multiple synchronization tongues are disposed 
along the lower boundary of the main synchronization 
tongue, i.e. their structure differs from the case shown in 
Fig. [1] (b) in the area of great values of coupling param- 
eter. 

Phase portraits of the first and the second oscillators 
computed inside the basic synchronization tongues are 
given in Fig. [4] Regime of doubled period on the basis of 
synchronous regime 2 is also shown in this figure. 

Non-identity of oscillators in the case of potential of 
the fourth degree may be introduced in one more way, 



namely by means of mutual displacement of the potential 
wells of interacting oscillators: 

CM*) = Y> U 2 (y)= { -V±^. (10) 

Parameter 6 controls in this case the position of the 
minimum of the fourth degree for the second oscillator. 
The dynamic regime chart is shown for this system in 
Fig. [5] We also observe the system of synchronization 
tongues. An interesting feature of synchronization pic- 
ture consists in this case in the fact that with the excess 
above the value S m 1 synchronization with the period 1 
displaces all other regimes and occurs even if the coupling 
parameter is indefinitely small. Phase plane portrait 
given in Fig. [2] indicates the fact that the first oscilla- 
tor dominates over the second oscillator (see scales in 
coordinate axis). 



V. CONCLUSION 

The analysis of synchronization in auto-oscillating sys- 
tems by means of the dynamic regime chart construction 
reveals an essential role of the anisochronism. If we in- 
sert nonlinearity of Duffing oscillator type, synchroniza- 
tion tongues become essentially wider and the situation 
of their overlapping occurs. Islands of doubled periods, 
" crossroad area" structures and chaotic areas appear in- 
side the synchronization tongues. For oscillators with the 
potential of the fourth degree we have observed the sys- 
tem of synchronization tongues, which depends on the 
parameter that controls relative steepness of the poten- 
tial wells of oscillators. The oscillation death effect dis- 
appears in this case and synchronization tongues modify 
themselves in the area of great values of coupling parame- 
ter. If mutual position of the potential wells of the fourth 
degree is changing, there is a critical value when synchro- 
nization with the period 1 occurs even if the coupling 
parameter is indefinitely small. 
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